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Abstract 

We introduce a method for extracting the cusp anomalous dimension at L loops from four-gluon 
amplitudes in J\f = 4 Yang- Mills without evaluating any integrals that depend on the kinematical in- 
variants. We show that the anomalous dimension only receives contributions from the obstructions 



introduced in hep-th/0601031 We illustrate this method by extracting the two- and three-loop 
anomalous dimensions analytically and the four-loop one numerically. The four-loop result was 
recently guessed to be / (4) = - (4Cf + 24( 2 C4 + 50C 6 - 4(1 + r)Cf) with r = -2 using integrability 



and string theory arguments in hep-th/0610251 Simultaneously, was computed numerically 



m 



hep-th/0610248 from the four-loop amplitude obtaining, with best precision at the symmetric 
point s = t, r = —2.028(36). Our computation is manifestly s/t independent and improves the pre- 
cision to r = —2.00002(3), providing strong evidence in favor of the conjecture. The improvement 
is possible due to a large reduction in the number of contributing terms, as well as a reduction in 
the number of integration variables in each term. 
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I. INTRODUCTION 



M = 4 super Yang-Mills (MSYM) is a very fascinating theory. Recently, important 
progress has been achieved in two different but connected aspects of the theory. On the one 
hand, insights from integrability have led to a proposal [l, 2| for all- loop quantities such as 



the cusp anomalous dimension f(g), which is the coefficient of the 



ogarithmic correction 



to the anomalous dimension of large spin twist-2 operators 3), 0, 0, |(3]. On the other hand, 
in [3] the integrand of the planar four-loop four-gluon scattering amplitude was found, and 
the amplitude was evaluated analytically to order 1/e 4 and numerically to order 1/e 2 in the 
dimensional regularization parameter e = 2 — D/2. 

The connection between these two remarkable results is that the coefficient of g 2L 
in the perturbative expansion of f(q) controls part of the IR singular behavior of the L- 



loop scattering amplitude 



lOt ] . Hence one can in principle extract that coefficient, 



and thereby determine f(g) order by order, by calculating four-gluon scattering amplitudes. 
However, in practice this is a very hard task. One reason for this is that the IR singularity for 
an L-loop amplitude begins at order l/e 2L , while the contribution from f(g) only appears in 
the 1/e 2 singularity. Therefore, the larger L, the higher one has to climb into the e expansion 
before uncovering the desired contribution from f(g). 

One of the results of |2j is a formula for f(g) whose expansion can be written in the form 

f(g) = Ag 2 - 4C 2 <? 4 + (4C 2 2 + 12&) g 6 - (4C 2 3 + 24^4 + 50C 6 - 4(1 + r)C 3 2 ) g 8 + 0(g 10 ), (1) 

where r was conjectured to have the value r = —2 and g 2 = gY M N/(87r 2 ). This formula 
refines an earlier guess in ll| which had the same form but r = 0. The non-zero value of r 
is due to a nontrivial phase factor in the spin chain S-matrix. Remarkably, the value r = — 2 



at four loops follows from a simple conjecture for the all-loop phase factor a. 
which has been shown to agree well with the AdS/CFT correspondence 12 



so given m 



13| 



On the other hand, the numerical evaluation of the 1/e 2 coefficient of the four-loop 
amplitude was carried out in The numerical calculation is done at different kinematical 
points x = t/s with results (r,x) given by (-2.028(36), 1), (-2.059(36), 2), (-2.062(45), 3), 
(—2.074(104), 15), where the quantities in parentheses denote the uncertainty in the last 
digits. 

An unfortunate complication inherent in this approach is that the L-loop cusp anomalous 
dimension is just a number, but the 1/e 2 coefficient of the L-loop scattering amplitude is, 



in general, a very complicated function of the kinematic invariant x = t/s. Isolating this 
single number is like picking a needle out of a haystack. For example ; the evaluation of the 



four-loop amplitude using Mellin-Barnes representations Il5l . Il6l . Il7| (the current state 
of the art) expresses the 1/e 2 coefficient as a sum of over 50,000 x-dependent terms 1 . Each 
of these terms is a multiple integral which in general is evaluated numerically. 

The fact that an appropriate combination of these terms conspires to add up to an 



x-independent number highlights an inefficiency in this calcu 



manifest in the fact that the numerical precision obtained by [7] depends on the value of x 
where the evaluation is performed, as noted above, even though the result of the calculation 
cannot depend on x. In terms of r, the precision reported by [7| varies from approximately 
1.5% at x = 1 to 5% at x = 15. 

In this paper we demonstrate that the cusp anomalous dimension may be calculated in a 
manifestly x-independent manner, and with a significantly smaller set of terms with reduced 
number of integrations variables. The starting point at four loops is the formula 

M (4) (x,e) = - (M (1) (x,e)) 4 - (M«(x, e))' M^(x, e) + 

+~ (M( 2 )(x,e)) 2 + i/WM«(Me) + 0(l/e), (2) 

which relates the four-loop amplitude to a polynomial in lower-loop amplitudes, together 
with a 'correction' proportional to f^ A \ Rather than calculatingboth sides of this equation 
as functions of x, and then reading off f^ 4 \ as was done in [7], we show that (and 
more generally f(g)) only receives contributions from the 'obstructions' discovered in flo| . 
Moreover, only the lowest order obstructions contribute. It turns out that obstructions can 
be systematically extracted from x-dependent integrals and the rest can be thrown away 
term-by-term, eliminating the need to evaluate all of them and then rely on non-trivial 



ational approach. This is 



cancellations. We carry out the relevant calculation analytically 



or the one-, two- and 



as a starting point 



three-loop amplitudes. At four loops we use Czakon's MB program 
and then implement the systematic extraction of obstructions. Quite nicely, it turns out 
that a large number of the integrals that come out can be reduced by simple applications 
of Barnes lemmas and corollaries. This process was also automated. After performing the 



1 Each "term" here refers to an integral of a single rational expression of T functions and their derivatives. 
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remaining integrations numerically, we obtain 

/( 4 ) = -117.1789(2), or r = -2.00002(3), 



(3) 



which is consistent with the prediction r = —2 within the precision of about 0.001%. 

In fact this computation of the cusp anomalous dimension is a simple application of a 
more general analysis of the structure of multi-loop amplitudes we present in this paper. 
Instead of working with amplitudes as functions of x, we work with their Mellin transforms to 
a new variable y. In y space, any amplitude can be uniquely expressed as a sum of two kinds 
of terms. The first kind are singularities, which take the form of ^-functions or derivatives of 
5-functions at y — 0. These singularities are the obstructions introduced in 19j; in x-space 
they correspond to polynomials in log 2 x. The remaining part of the amplitude is smooth 
at y — 0. We show that the obstructions in L-loop amplitudes must separately satisfy any 
polynomial relation of the form (T5]) and that furthermore the cusp anomalous dimension 
only receives contributions from the leading singularity 5(y) at order 1/e 2 . To summarize: 

Our prescription for computing the cusp anomalous dimension is to read off the 
numerical coefficient in front of S(y)/e 2 on both sides of the Mellin transform of 
equation (T5]). 

This paper is organized as follows. In section 2 we review how the cusp anomalous 
dimension may be extracted from the infrared singularities of gluon scattering amplitudes. 
In section 3 we define obstructions in detail, discuss their most important properties, and 
explain how they may be used to calculate the cusp anomalous dimension. In section 4 we 
illustrate the method at the two- and three-loop level exhibiting the form of the obstructions 
explicitly. In section 5 we give some details on the computation of the four-loop cusp 
anomalous dimension from obstructions. In section 6 we summarize our conclusions and 
list some possible promising directions for future work. Finally in appendix A the reader 
familiar with Mellin-Barnes representations will find detailed instructions for using Czakon's 
MB program to extract the obstructions from an amplitude. 
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II. STRUCTURE OF THE L-LOOP FOUR-GLUON AMPLITUDE 



We let A^ L \e, s, t) denote the planar L-loop four-gluon amplitude in M = 4 super- Yang- 
Mills theory, dimensionally regulated to D = 4 — 2e dimensions. We further define 



M< L \x,e) = {sty L / 2 A^ L \s,t,e)/A^\s,t), x = t/s. 



(4) 



It is convenient to divide by the tree-level amplitude because supersymmetry determines the 
helicity structure of the L-loop amplitude to be the same as that of the tree-level amplitude; 
therefore the ratio is a function only of e and the Mandelstam invariants. The prefactor 
(st) eL l 2 is convenient because the resulting quantity is a function of the ratio x = t/s 
only, and satisfies the symmetry 



M (i) (x,e) = M (L) (l/x,e). 



(5) 



However we caution the reader that it is not standard in the literature to factor out (st) eL//2 
as we have done. 



Expressions for M^ L \e, x) in terms of a basis of integrals are known for L = 1, 2, 3, 4 ( 20], 



21 1, 0] respectively). For L < 4 the basis is determined by the structure of the simplest 
unitarity cuts, i.e, double cuts. This analysis led to the "rung rule" which basically states 
that the L-loop amplitude can be constructed from the diagrams that give rise to the (L — 1)- 
loop amplitude by adding one rung to all diagrams in such a way that no internal bubbles 
or triangles are produced 21]. For L > 4, the rung rule gives only the part of the amplitude 
that can be determined by double cuts. The four-loop amplitude is then special in the sense 
that it is the first order in which the rung rule is incomplete. 

In a remarkable recent effort [7J the two diagrams missing in the four-loop amplitude were 
found. There it was noticed that these two diagrams, as well as the ones produced by the 
rung rule, satisfy some intriguing conformality properties [22j, |23j that have been studied 



developed in [25 



recently in 24 L This new insight, together with some generalized unitarity arguments 



27 



281 ] . might lead to a systematic way of generating the diagrams at 
any order in perturbation theory. Very recently a proposal for the integrand of the L = 5 
loop amplitude was presented [29]. In Figure 1 we show the integrals contributing to the 
L < 4 loop amplitudes. 
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FIG. 1: Here we depict the topologies of the integrals which contribute to the L-loop four-gluon 
amplitude, for L < 4. The numerical coefficients and the numerator factors in the integrand are 
given by the rung rule for all but the last two four-loop diagrams. We refer the reader to {7] for 
further details. 



A. Infrared Singularity Structure 



Let us now review the known IR singularity structure of scattering amplitudes of gluons 
in Af = 4 SYM. The cusp anomalous dimension controls part of the 1/e 2 singularity and 
therefore it can be computed once the amplitudes are known explicitly. 

In particular, it is known 0, Q, Q that the infrared singularities of an L-loop gluon 
amplitude, which contain l/e 2i and lower poles, can be expressed as a polynomial in lower- 
loop amplitudes according to the formula 



M( L \x,e)=XW [M<"(i,e)] + ^f {L) M^(x, Le) + 0(1/ e) 
where the polynomials X^ are conveniently encoded in the formula 



(6) 



X (L) [M w (z,e 
The one-loop amplitude takes the form 

M (1) (x,Le) 



M {L) (x, e) - log ( 1 + jr a l M (l) (x, e) 



(7) 



a L term 



L 2 t 2 



+ 0(1). 



(8) 



Therefore, (jSJ) implies that the L-loop cusp anomalous dimension is given by the following 
deceivingly simple formula 



= _ 2L 2 r M (i)( X;C ) _ X (L) r M (0( x?e ) 



1/e 2 



(9) 



where the subscript on the right instructs us to take the coefficient of the 1/e 2 singularity. 



Note that is a number; therefore the right hand side of must be x independent. 
The calculation of \\ does not make this x-independence manifest, so was computed for 
various values of x as a consistency check, as mentioned in the introduction. The cancellation 
of the x dependence in fl9]) is precisely the delicate step that we manage to avoid with our 
obstruction technique as we explain in the next section. 

Relations similar to ([6]) hold in any massless gauge theory, but in J\f = 4 Yang-Mills 
theory it is believed that the stronger statement 



M (i) (x, e) = X (L) [M {1) (x, e)] + ^f {L) (e)M« (x, Le) + C (L) + 0(e) 



(10) 



holds, where are numerical constants and /^(e) are functions of e only, with /^(0) = 
/( L ) being the L-loop cusp anomalous dimension. This relation was originally proven for the 
two- loop four-particle amplitude in 30J], and has been subsequently shown to hold for the 



31 



32|. 



three-loop four-particle amplitude [18( as well as the two-loop five-particle amplitude 

Below we will make use of ffTU]) at two and three loops, as well as <Q at four loops, so we 
record here the explicit forms of these relations, 



M^(x,e) 
M (3) (x,e) 
M (4) (x,e) 



(x,2e) + C (2) + 0(e), 



-(M«( 



+ M (1) (e, x)M (2) (x, e) + ^/ (3) (e)M (1) (x, 3e) + C (3) + 0( 
j (M«(x,e)) 4 - (M«(x,e)) 2 M®(x,e) + M^(x, e)M^(x, e) 
+\ (M( 2 )(x,e)) 2 + V)M«(x,4e) + 0(l/e), 



where pj, [3fl 



(11) 



r/ 



(2)/ 



C (2) 



C(3) 



-C 2 -C 3 e-C4e 2 + 0(e 3 ), 

yC4 + e(6C 5 + 5C2C3) + 6 2 ( Cl C 6 + c 2 C 3 2 ) + 0(e 3 ) 
341 2 \ . / 17 2 \ 2 



216 9 



9 9 



(12) 



The quantities c\ and c 2 are currently unknown, but they drop out of ffTUl) . The quantities 



/^(0) agree with the original 



loops carried out in 



33 



34 



35|. 



calculations of the cusp anomalous dimension up to three 
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III. OBSTRUCTIONS AND THEIR APPLICATIONS 



In this section we first refine the definition of obstructions given in 19J, clarifying their 
more formal mathematical structure. This allows us to give an all order in e formula for the 
one-loop obstruction. We then explain the systematic way of implementing the extraction of 
obstructions at higher loops using as a starting point the program MB by Czakon 17| • Finally, 
we give the formula for the L-loop cusp anomalous dimension in terms of the obstructions 
at < L loops. 



A. Obstructions — Theory 

The state-of-the-art technique for calculating L-loop integrals makes use of Mellin-Barnes 



representations Q, [l5| (see [l6| for a thorough introduction). This technique leads us to 
work not with an amplitude (x, e) but rather with its Mellin transform (y, e). These 
are related by 

M^(x,e) = P°° £L x v F<%,e). (13) 

It is fairly simple to write down an expression for F( L \y,e) using the Feynman rules for 
any desired L-loop diagram. However these expressions are unwieldy because they typically 
involve multi-dimensional integrals of ratios of gamma functions. For example, the eight 
diagrams which contribute to the four-loop amplitude require Mellin-Barnes representations 
of dimensionalities from 10 to 14. 

In the representation (fTBl one typically cannot expand around e = under the integral 
sign because the function F has singularities which collide to pinch the integration contour 
as e — > 0. We isolate these singularities using the formula 

lim — !— = V-± tt%) (14) 
w-*o y ±lu y 

and the generalization to higher singularities, given by taking derivatives of this equation 
with respect to y. The formula is missing some familiar factors of i because the variable y 
runs along the imaginary axis, and hence we take u real 2 . 



2 It is subtle to make sense of (^-functions of complex arguments, and we make no attempt to do so here. 
Rather, we emphasize that in the analysis of this section we keep the y contour precisely on the imaginary 
axis, meaning that one could rewrite all of the formulas in this section in terms of the real variable z = —iy 
if one preferred. In this case one would have to take w imaginary in |T] 



S 



Isolating all of the singularities at y = leads to a formula of the form 

M (L) (x,e) = l +t0 ° x y H^ L \y,e) + P {L \\ogx } e), (15) 
J-ioo 2m 

where -f denotes a principal value integral, and we can now expand in e under the integral 
sign. Since the Mellin transform of -^5(y) is proportional to log fc a;, the residual terms 
called in (P) are guaranteed to be polynomial in log a; (order by order in e). These 
terms were called 'obstructions' in [l9] because their presence indicates a failure to be able to 
assemble everything under a common y integral at e = 0. We see therefore that obstructions 
correspond to 5(y) -function singularities in the Mellin transform of an amplitude. In fact, 
the symmetry ([5]) guarantees that P^ L ' (log x, e) will be an even polynomial, and it also 
guarantees that H^ L \y, e) is an even function of y. 

Let us illustrate this idea by calculating the obstruction to all orders in e for the one-loop 
amplitude. This result is useful since the 1 /e 2 term of (x, e) is controlled in part by the 

£ 2L-4 term of M (l) 

(x, e) as can be easily seen from ([7]). 
In terms of the function 

F {1) (y, e) = -^pT^yr(l + f e + y)T\-\e + y)T\-\e - y)T{\ + f e - y), (16) 

the one-loop amplitude is given by 

M (1) (x,e) = ^L x v F^(y,e) } -2 < e < 0. (17) 

J-ioo 2m 

As e — > 0, the two poles at y = ±|e collide with the integration contour. We can put those 
poles directly on top of the integration contour, with a principal value prescription, at the 
expense of subtracting off one-half of the appropriate residues as in ( 1251) . Therefore we find 

Art, j I 

JLx»FM(y,e) + -Res y=+€/2 [xyF^(y,e)} - -Res^/^F^y, e)]. 
-too Zm Z z 

(18) 

The quantity F^(y,e) under the integral is to be expanded in e inside the integral. The 
relative sign between the last two terms arises because the pole at y = +e/2 approaches the 
integration contour from the left while the y = — e/2 pole approaches the contour from the 
right. The obstruction is the sum of these two residues. Using the explicit formula ffTBT) . we 
find that the one-loop obstruction to all orders in e is 

P (1) (logx, e) = ^ [{x +t/2 + x" e/2 )(^(l + e) - 2^(-e) + ^(1)) + (x +t/2 - oT e/2 ) log re] , 

(19) 




(20) 



As expected based on our general arguments, the e expansion of the obstruction fTTTJT) is 
an even polynomial in logx. We have also obtained an analytic formula for the two- loop 
obstruction to all orders in e, but the expression is quite lengthy and involves some remaining 
Mellin-Barnes type integrals that we do not know how to perform analytically at finite e. 

B. Obstructions — Implementation 

In the previous subsection we gave a formal definition of obstructions and used it to 
compute the one-loop obstruction to all orders in e. However, at higher loops the calculation 
rapidly becomes overwhelming, and must be implemented in a computer algebra system. 
Fortunately we can use Czakon's MB program which is designed to aid in the manipulation 
of Mellin-Barnes integrals as the starting point of the implementation. Let us now outline 
the steps we took in order to use MB to calculate obstructions. This outline will probably 
only be of interest to those with some familiarity with Mellin-Barnes integrals, or with MB 
in particular. A more concise outline of this procedure is given in the appendix. 

We begin with a Mellin-Barnes representation for each integral / contributing to the 
desired amplitude, 



where f(zi,e) is a linear function of the z^ and e, and F{z^e) is in general a complicated 
ratio of T functions whose arguments are linear in the Z{ and e (a four-loop example is given 
below). The integration variables run from — ioo to +200 along contours that can be chosen 
to be parallel to the imaginary axis, with appropriate real parts. 

The representation (I2T!) of I(x, e) generically is not valid in a neighborhood of e = 0, 
rather one has to start with a representation valid around some e = eo < 0. One then 
analytically continues the formula to a neighborhood of e = 0, where one is free to develop 
the Laurent expansion of the amplitude that will exhibit all of the IR poles. In the process 
of this analytic continuation one must cross poles of the various T-functions in F(zi, e). The 
end result is the sum of the original integral, now valid in a neighborhood of e = 0, plus a 
sum of residue terms which have fewer integration variables than the original integral. 




(21) 
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Once this process is completed, the coefficient at any desired order in e is given as a sum 
of multiple integrals, most of them containing a factor of x^ Zi ^ to a power that is linear in 
the remaining integration variables. Some other integrals do not depend on x, or their x 
dependence is an overall factor which is a power of log x. 

In each of the integrals which contains a factor of x^ Zi ' we perform a linear change of 
variables f{zi) —> y. The next step is to take all integrals that contain an x y factor and 
combine them under a single y integral sign. Generically, this is impossible because the 
integration contour for y can be different in each individual term. For example, two typical 
contributions might look like 



/•/Ji+ioo Ay r/3 2 +ioo Ay 

h(x) = / ^ xTMHxiy), h(x) = / £L xyT(y)H 2 (y) (22) 



where f3\ > > (3 2 are on different sides of y — 0. Therefore there is an obstruction at y = 
that prevents us from combining these two integrals. The obstruction is a simple pole in 
this case. Therefore, we can write 



r[32+tco rty 

h(x)+I 2 (x)= / -JL x y T (y)(H 1 (y) + H 2 (y)) + A (23) 

where A is the obstruction, i.e, the contribution obtained from the residue at y = 0, 

A = ^-i xyT{y)H x {y). (24) 
2m J y =o 

As mentioned before, the presence of these obstructions that prevent one from naively 
combining all terms under the same integral sign were discovered in [3]. In that paper, 
the goal was to prove relations of the form (TTTJT) as an identity under the integral sign. 
Therefore, obstructions were undesired features of the Mellin-Barnes representation and 
differential operators were designed to annihilate them. Since our main goal in this paper 
is to compute the cusp anomalous dimension, it is precisely the obstructions that we are 
interested in. 

In 19), we were content to collect all terms under a single with an arbitrary contour (3. 



However, as we have already indicated in the previous subsection, it turns out that the most 
natural and useful notion is to take (3 = so that the y contour sits along the imaginary axis, 
with a principal value prescription. If (3 is infinitesimally small — which in practice means 
that it is sufficiently small to ensure that there are no more poles between the y contour and 



11 



the imaginary axis — then the final step of pushing the y contour to the imaginary axis can 
be accomplished via 

J/3-ioo 2lXl J-ioo 2m 2 

The final step in calculating the full obstruction in any amplitude M(e, x) is to add up 
all of the contributing integrals J(e, x). The above steps provide a constructive proof that 
any amplitude (and moreover, each individual integral I(x, e) contributing to an amplitude) 
can be written in the form 

«(y,e)+P(logx,e) (26) 

-ioo 2711 

for some H(y,e) and P(logx, e). Note that the sum over contributing integrals I(x,e) 
imposes that M(e, x) satisfies the symmetry ([5]). This ensures that the resulting H(y, e), after 
all contributions are added together, satisfies H(y,e) = H(—y,e), and that the obstruction 
P(logx, e) must be symmetric under x <-> l/x, so it is an even polynomial in logx. 

C. Cusp Anomalous Dimension Prom Obstructions 

Finally we explain how to exploit obstructions to simplify the calculation of the cusp 
anomalous dimension using equation ([9l). There are two essential properties of obstructions 
that enable this calculation. 

The first essential property is that obstructions satisfy a product algebra — i.e. if Pi and 
P2 are the obstructions in two amplitudes Mi and M 2 , then the obstruction in the product 
M\(x, e)M 2 (x, e) is simply the product of obstructions Pi (logx, e)P 2 (logx, e). 

This is easy to see in y-sp&ce, where multiplication of Mi(x,e) becomes convolution of 
the corresponding Fi(y, e). Delta functions (or derivatives of delta functions) convolved with 
other delta functions still gives delta functions. However, delta functions convolved with a 
smooth function give a smooth function, as does of course a smooth function convolved with 
a smooth function. 

This fact allows us to calculate, order by order, the obstructions in the 1,2,3, etc. loop 
amplitudes once and for all, and then they can subsequently be plugged into ([9]) to study the 
cusp anomalous dimension. Further more, any polynomial iteration relation of the form (flOl) 
that is satisfied by multi-loop amplitudes must also be satisfied by their corresponding 
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obstructions. This is just the statement that we can take the Mellin transform of both sides 
of the equation and truncate to the terms proportional to (derivatives of) delta functions at 
y = 0. 

The second essential property is that the cusp anomalous dimension only receives con- 
tributions from obstructions, and in fact only from the leading obstruction proportional to 
8(y) (as opposed to -^5(y) for k > 0). To prove this property, let us recall formula (jHJ) for 
the L-loop cusp anomalous dimension 



f(L) = _ 2L 2 [ M (L) (x,e) -X (L) [M (/) (x,e) 
Use the decomposition (1261) to write 



l/e 2 



(27) 



M (i) (x,e)= / °° ^ x y v\H {L) (y,e)] + P (i) (logx, e) 
J-ioo 2m L J 



(28) 



and 



X™ [M« ( X , e)] = / l °° £L x*V \K& (y, e)l + Q (L) (1 ogx, e). (29) 
for some K and Q. In fact the product algebra structure tells us that Q^ L \logx, e) is given 



in terms of the obstructions in lower loop amplitudes by 



Q (L) (log x,e)=XM [P {1) (log x, e)] . 



(30) 



Here we write the principal value integral as the integral of the principal part, with a 
conventional integral along the imaginary y-axis instead of The reason is that we would 
like to collect all of the terms in equation (f2"T)) under the single y integral 



dy 



^L x *5(y)fW = -2L 2 [ l °° JLxyv\HM(y,e)-K( L \y,e) 
too 2m J-ioo 2m 1 



-21} 



IOO 

ioo 



1/e 2 



dy 
2m 



x y 



Q 



d_ 

dy' 1 



8{y). (31) 

l/e 2 



Note that we use 5-functions that contain an extra normalization factor of 2m' to compensate 
for our normalization of the y integral. The e expansion is performed inside the y integral, 
as explained in the previous section, and the l/e 2 coefficient of the term on the second line 
is a polynomial in d/dy acting on 5{y) inside the integral. 

Because equation (13TT) must hold for all x, the integrand must satisfy the equation point- 
wise in y. This means that the principal parts and the 5-function parts must separately 
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satisfy the equation, and in particular the coefficient of ^5{y) must satisfy fl3Tl) separately 
for each k. By reading off the coefficient of S(y) on both sides of (!3T1) we conclude that 



./ 



-2V 



(0,e)-Q^ (0,e) 



1/e 2 



(32) 



Finally, using fl30l) we arrive at our main formula for computing the L-loop cusp anomalous 
dimension from the obstructions of I < L loop amplitudes, 



f(L) = _ 2L 2 \ p {L) ( , e ) -X (L) [P w (0,e)] 



l/ £ 2 



(33) 



IV. ANALYTIC RESULTS AT ONE, TWO AND THREE LOOPS 



In this section we illustrate our method by showing that an analytic computation of 
and is possible without computing any x-dependent integrals. In order words, 
none of the technology available for evaluating Mellin-Barnes integrals in terms of harmonic 
polylogarithms is needed. We also show that relations like (fTOl) that are expected to hold 
for the full amplitude must and do hold for the obstructions alone. 

A careful analysis of the one-loop amplitude was carried out in section III. A and the 
formula (fT9l) displays the obstructions in this amplitude to all orders in e. The computation 
of only requires the expansion of P^(logx, e) through order e 2L ~ 4 . Since in this paper 
we are interested in L < 4 it is sufficient to expand P^ through order e 4 : 

"irCs" 



pW(i 



ogx,e) 
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h — — 
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70CI 



+ C(e 5 



(34) 



4320 9 5760 4608 

We have also obtained an exact formula for the obstructions in the two-loop amplitude. 
However instead of displaying the complete complicated formula we show its expansion 
through order e 2 , which is sufficient for the purpose of isolating its contribution to f^: 

5tt 2 
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ogx, e 



+ 



2 1 
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At three loops we have worked out the obstructions through finite order, finding 
P (3) (logx,e) = 



1 
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1617T 4 
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1632960 18 576 256 256 

It is a simple but delightful exercise to check that these three expressions obey the two- and 
three-loop iterative relations (1101) . (1121) ■ Although amusing, it is of course guaranteed from 
our general discussion of the properties of obstructions in the previous section. 



V. FOUR-LOOP CALCULATION 



In this section we explain the details of our numerical computation of the four-loop cusp 
anomalous dimension using obstructions. 

The first step is to find an expression for the four-loop amplitude in terms of a small 
number of integrals. As mentioned in section II, the rung rule does not give the full amplitude 
for L > 3. Finding the remaining contributions turns out to be a very hard problem. Luckily, 
in [?| , the new non-rung-rule contributions were found and the full answer passed nontrivial 
consistency checks like unitarity and the correct IR singular behavior. Here we take the 
diagrammatic formula of [7|, shown in Figure 1, as our starting point. 

The second step is to find Mellin-Barnes representations for each of the integrals con- 
tributing to the amplitude. This is a straightforward process but such representations are 
not unique. Here we give as an illustration a 10-fold Mellin-Barnes representation for the 
four-loop ladder integral (the tenth integral / x y is suppressed) 



F(y,e) 



3 4<E7 



9 dz- 

II ^ T(y - 2e)r(l + 2e - y)T(-l - e - Zl )Y{- Zl )Y{-z 2 ) 



r(-2e) - l=1 

xr(l + z x + z 2 )Y(-l - e - z 3 )T(2 + e + z x + z^)Y{l - z 2 + z 3 
xr(l - Zl -z 2 + z 4 )r(-2 - 2e - z 5 )r( Zl - z 5 )T(-z 4 + z 5 ) 
xT(— 3 - 2e + z 2 - z 3 - z A - z 6 )T(-z 6 )T(3 + e + z 3 + z 5 + z e ) 
xr(-3 -2e-y- z 7 )T{l -z A + z 7 )T(-3 - 3e - z 8 )T{z 5 - z 8 ) 
xT(4 + 2e + y + z$)T(-z 7 + z$)T(-l - e + z 4 + z 6 - z 7 - z 9 ) 
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xr(-3 -2e-y-zg- z 9 )T(-z 9 )T(3 + e + y + z 7 + z 9 ) 
xT(l - z 5 - 2 6 + Z 8 + Zq)/ 

(r(-l - 3e - z 1 )F(l - 2 2 )r(3 + e + z x + z 3 )Y{-2 - 4e - z 5 ) 

xr(i - z 4 + * 5 )r(i - z 6 )r(-3 - 5e - z 8 )r(i - z 7 + * 8 )r(i - z 9 )). 



(37) 



The reason we display this one is that this representation has one less integration variable 
than the one given in , reflecting the non- uniqueness of Mellin-Barnes representations. 

The third step is to extract the obstructions at order 1/e 2 from the Mellin-Barnes inte- 
grals. This step gives rise to a list with a large number of integrals. All these integrals are 
independent of x. 

The fourth step is the simplification process. It turns out that a large fraction of these 
obstruction integrals can be reduced to integrals with a smaller number of integration vari- 
ables by applications of Barnes lemmas and their corollaries (see appendix D of for a 
nice collection of lemmas). 

The result of this process is a formula for in terms of 13 5-fold integrals, 131 4-folds 
and many more 3-, 2- and 1-folds. As an illustration we give the following "analytic" formula 
for the four-loop cusp anomalous dimension: 
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+ 



57661 r dz x 
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(38) 



In this expression the integration contours for each variable are straight lines parallel to the 
imaginary axis with real parts given respectively by 
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(39) 



The final step is the numerical evaluation of these integrals, which we performed using 



CUBA's Cuhre algorithm 36] as implemented by Czakon's MB program. It is important to 
mention that the main source of numerical error comes from 5 and 4 folds. The full list of 
integrals is available from the authors upon request. 

In the simplification process we implemented only the first and second Barnes lemmas 
and some of their simplest corollaries. We believe that a more detailed analysis will show 
that some major simplifications are still possible and a fully analytic answer is within reach. 
We believe that this is a very important problem but we leave it for future work. 



VI. CONCLUSIONS AND FUTURE DIRECTIONS 



In this paper we clarified the nature of 'obstructions' in L-loop four-gluon amplitudes 
in M = 4 Yang-Mills theory. Their name comes from the fact that they are the terms 
which obstruct writing an amplitude as an inverse Mellin transform of a smooth function 
near e = 0. These obstructions were shown to possess all the information needed for the 
computation of the cusp anomalous dimension at L-loop order. 

We implemented a systematic way of extracting obstructions, thus leading to a new 
method for computing the cusp anomalous dimension. As an application of our method we 
recovered the known two- and three-loop results already given in the literature 
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30 



34 



35) 



33|, 



We also used the recently obtained [/J integrand for the four-loop four-gluon amplitude, 
given by the diagrams shown in Figure 1 , to compute numerically the four- loop cusp anoma- 
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lous dimension. Our result is 

/( 4 ) = -117.1789(2), (40) 

where the quantity in parentheses denotes the uncertainty in the last digit (as reported by 
the CUBA numerical integrator). 

It would be highly desirable to find an analytic formula for /W. This is something which 
is definitely within reach and we leave it for future work. However, we can try to see how far 
one can get with (1401) and some additional assumptions. The first assumption one can make 



is the KLOV transcendentality principle 35|, |37|, which assigns a degree of trascendentality 



k to (k, 1 to it and zero to rational numbers. The degree of a monomial is computed as the 
sum of the degree of each factor. The transcendentality principle then states that is 
homogeneous of degree 2L — 2. This means that must be a linear combination of 7r 6 , £f , 
and 7r 3 (s with rational coefficients. Clearly, if we take this principle literally there is nothing 
we can conclude from (l4"0~j) for there would be an infinite number of possibilities consistent 
with it. 



Let us take a stronger form of the principle which was observed in integrability studies [11 ]. 
The idea is to allow only numbers to appear in the monomials, and only with integer 
coefficients with sign (— 1) L+1 . Then at four loops we would have the ansatz 

/ (4) = -KC 2 3 + ™ 2 C 2 C4 + n 3 ( 6 + n 4 Q) (41) 

with four positive integers n,. There is of course ambiguity in the choice of (ni,n,2, n^) since 
7^6 = 4C2C4 an d 5C2C4 = 2C|- Then we can ask how many inequivalent choices (ni, n 2 , n 3 , ri 4 ) 
are consistent with our result (140]) within error bars. It turns out that there is only one 
choice, namely 

- (4C 2 3 + 24C2C4 + 50Ce + 4C|) ~ -117.1788285, (42) 
which is the value conjectured in j^j]. The next nearest value 

- (2C 2 3 + C2C4 + Ce + 73Cf) « -117.1801235 (43) 

is over six sigma away from our numerical data. For comparison, we note that even with the 
best precision obtained in p| there are more than 200 inequivalent choices of (n\, n 2 , n 3 , n^) 
within the error bars. 

Another interesting problem for the future is to study the feasibility of a higher loop 
calculation. As suggested in jl], if the integrability result holds to all orders in perturbation 
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theory, then there must be some hidden structure in field theory that allows the same power 
of computation. We hope that the theory of obstructions elaborated in this paper can 
be refined even more so as to produce a technique applicable in practice to all orders in 
perturbation theory. There are some indications to be hopeful about this. One is the fact 
that we were able to obtain a closed formula for the one loop obstruction to all orders in e 
(see [38] for similar all-orders quantities). The second is that, assuming the transcendentality 
principle, the anomalous dimension should only receive contribution from the highest degree 
terms. The rest must cancel out. It turns out that in our calculation many integrals gave 
rise to terms of the "wrong" degree which later on cancel. It is conceivable that the terms of 
wrong degree might be isolated and thrown away systematically leaving a simple expression 
for the highest degree terms — mirroring the way in which our method of obstructions isolates 
and systematically throws away the x-dependent terms which we also know must cancel. We 
also leave this fascinating question for the future. 

It is also very tempting to suggest that a five-loop calculation might be feasible with our 
technique or a further refinement of it. Something certain is that a five- loop computation 
of the anomalous dimension using the IR structure of the five-loop amplitude 29J] and 
then directly using MB to numerically evaluate both sides of (j2J) is well beyond our current 
computational power (the numerical computation of the four-loop anomalous dimension 
presented here already required over a CPU- month on 3 GHz Xeon computers). 

It would be interesting to understand the connection between obstructions and the con- 
verse mapping techniques employed for example in [39] to find asymptotic expansions of 
amplitudes, and to see whether any of the simplifications we obtained can be applied to the 
direct calculation of the cusp anomalous dimension from three-point functions (as in [dpi]). 

Finally it would be interesting to investigate the applicability of similar techniques to 
other theories, including possibly QCD. 
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APPENDIX A: IMPLEMENTATION DETAILS 

Here we provide a concise outline explaining a straightforward way to harness the power 
of Czakon's MB program to calculate the L-loop cusp anomalous dimension using the method 
of obstructions as outlined in the main text of this paper. Recall that to calculate the cusp 
anomalous dimension we first need to extract the numerical coefficient of 5(y)/e 2 in the 
Mellin transform of the L-loop amplitude. This may be accomplished as follows: 

1. For each diagram that contributes to the L-loop amplitude, find a valid Mellin-Barnes 
representation and run MB to read off the coefficient of the 1/e 2 pole. The result is a 
list X of x-dependent integrals valid at e = 0. In each resulting term the x-dependence 
can take four possible forms: (a) independent of x, (b) proportional to a power of 
log x, (c) x to the power of a linear function of the remaining integration variables, or 
(d) as in (c) but with a prefactor of some power of log x. 

2. Throw away any integrals in X that contain any power of logx — these terms cannot 
contribute to the cusp anomalous dimension. 

3. In the terms containing x to some power, make a linear change of variables to a new 
variable y to write the x dependence in each term as x y . Of course, the location of 
the y contour will in general be different in each term. 

4. Take each integral which contains x y and perform an analytic continuation (by running 
the MB program again) to shift the y contour infinitesimally close to y = 0. Here 
'infinitesimally close' means, in practice, close enough to guarantee that there are 
no remaining poles between the y contour and the imaginary y-axis. This procedure 
produces a new list of ^-dependent integrals Y. The x-dependence in each term is either 
a power of logx, x y , or independent of x. 
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5. Throw away any integrals in Y that contain any power of log a?. 

6. Finally, for each term that still contains / dyx y , push the y contour onto the imaginary 
axis using (125]) and throw away the principal part as well as any terms containing logx. 
The result is a list of x- independent integrals Z. 

7. Add together the integrals Z obtained from all of the diagrams contributing to the L- 
loop amplitude. The resulting set of integrals adds up to a number . The interpreta- 
tion of this number is that it is the coefficient of ^S(y) in the Mellin transform of the L- 
loop amplitude M^ L \x,e). For example, at four loops we find ~ —174.285193(6). 

The L-loop cusp anomalous dimension is then given in terms of the number s^ L ' and the 
obstructions in the lower-loop amplitudes according to the formula fl33|) 



f L ) = -2L 2 [>>-XW [P<')(0,. 



(Al) 



where is the polynomial defined in ([6]) above. 

The full L-loop obstruction (including all of the log fc x terms) can be computed without 
too much additional work. One simply follows the above outline, order by order in e, but 
without throwing away any log fe x terms. 

The above procedure could actually be simplified somewhat by making an appropriate 
change of variables to y and choosing the contour for y to lie along the imaginary axis 
in the original Mellin-Barnes representation, before running MB. One would then use MB to 
push e infinitesimally close to e = 0, and then finally use formula ( |25|) for the final push to 
e = 0, obtaining a principal value integral plus obstructions. Note that this is how we chose 
to present the calculation of f|T9l . Although conceptually simpler, we did not employ this 
method at four loops because some experimentation revealed that this method seemed to 
give rise to a larger number of terms than the above procedure. 
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